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ABSTRACT 

Context. A possible correlation between CO luminosity and its line width (FWHM) has been suggested and denied in the 
literature. Such claims were often based on a small, or heterogeneous sample of galaxies, and thus inconclusive. 

Aims. We aim to prove or dis-prove the LJ-^-FWHM correlation. 

Methods. We compile a large sample of submm galaxies at z > 2 from the literature, and investigate the L'CO-FWHM relation. 
Results. After carefully evaluating the selection effects and uncertainties such as inclination and magnification via gravitational 
lensing, we show that there exist a weak but significant correlation between and FWHM. We also discuss a feasibility to measure 
the cosmological distance using the correlation. 
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1. Introduction 


A possible correlation between CO luminosity 
line width at high-z has been suggested in the literature. 
Both well et al.l (12013h observed 32 submm galaxies (SMG) at 
1.2 < z < 4.1 in '^CO and suggested a correlation between CO 
luminosity and line width (FWHM) in high-redshift starbursts, 
likely relating to baryon-dominated gas dynamics within the re¬ 
gions probed (See their Fi g. 5). Using 1 5 SM G at z=2-4 with 
known lens magnifications, Harris et ^ (l2012l) showed that the 
SMGs fall close to a power-law with small scatter in CO lumi¬ 
nosity and line width (FWHM) plane (See their Fig.7). 

Bright SMGs are often gravitationally-lensed. Such relation 
between and FW HM can be used to estimate magnifica¬ 
tion factor (p). Indeed, [Harris et'aP (1201 2h use d the relation to 
estimate luminosity of SMGs with unknown p. iFu et al.l (1201 3h 
also used this relation to estimate a small p (~1.8) of merging 
galaxies at z=2.3 to conclude their st ar formation rate were 2000 
Moyr^^ (See their Fig. S7). Similarlv. lRiechers et alJ(l2013l) used 
the relation to estimate p (1.5+0.7) of a starburst galaxy at 
z=6.34, finding its SFR of 2900 MQyr~^. 


Conversely, ICarilli & Waited (1201 3l) examined the relation 
using a much larger sample, and found no significant correla¬ 
tion between and FWHM (See their Fig.5). However, their 
sample was heterogeneous including various types of galaxies 
such as SMM, QSO, radio galaxies, 24jum sources and so on. 
Therefore, the existence of LJ-^-FWHM correlation is still con¬ 
troversial. 


In this paper, we compile a larger sample of z > 2 galax¬ 
ies with CO detection from the literature, to re-examine the 
Lpg-FWHM correlation at high-z. We carefully examine vari¬ 
ous sources of possible uncertainties such as different types of 
galaxies, luminosity conversion from different CO transitions. 



Fig. 1. Correlation between L’^glp and FWHM for various samples. 
The blue dotted line connects the median in equally spaced logarithmic 
bins. The black solid line is the best-fit power-law. 


and errors in estimating p. Unless otherwise stated, we adopt a 
cosmology with {h, Q.,„, Qa) = (0.7,0.3,0.7). 
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Table 1. Correlation between L'^qIp and FWHM for various samples. Dispersion is around the best-fit power-law. p-value of <0.05 can be 
considered significant. 


Sample 

N 

Spearman’s coeff. 

p-value 

Kendall’s coeff. 

p-value 

Dispersion (km/s) 

All 

195 

0.55 

1.7x10-*'’ 

0.39 

< 1.0 X 10-**’ 

187. 

Submm galaxies 

72 

0.53 

2.1x10-'’ 

0.37 

4.2x10-*’ 

272. 

QSO hosts 

31 

0.53 

2.0x10-3 

0.39 

1.9x10-3 

245. 

Submm galaxies -hQSO hosts 

103 

0.52 

1.4x10-3 

0.37 

< 1.0 X lO-*® 

218. 

CO(4-3) only 

26 

-0.34 

9.2x10-3 

-0.21 

1.4x10-* 

279. 

CO(3-2) only 

62 

0.12 

3.5x10-* 

0.09 

2.8x10-* 

217. 

CO(2-l)only 

38 

0.18 

2.8x10-* 

0.14 

2.2x10-* 

182. 

CO(l-O) only 

43 

0.22 

1.5x10-* 

0.15 

1.6x10-* 

121. 

All but CO(4-3) 

143 

0.52 

4.0x10-** 

0.37 

< 1.0 X 10-**’ 

206. 

All but CO(3-2) 

107 

0.56 

2.5x10-*** 

0.40 

< 1.0 X 10-*** 

201. 

All but CO(2-l) 

131 

0.58 

5.9x10-*3 

0.42 

< 1.0 X lO-*® 

231. 

All but CO(l-O) 

126 

0.52 

4.8x10-*'* 

0.37 

< 1.0 X lO-*® 

237. 

p known 

85 

0.42 

7.3x10-3 

0.30 

5.3x10-3 

156. 

p unknown 

no 

0.47 

2.4x10-3 

0.33 

3.6x10-3 

205. 

24pm sources 

8 

0.45 

2.6x10-* 

0.36 

2.2x10-* 

204. 

BzK galaxies 

4 

0.80 

2.0x10-* 

0.67 

1.7x10-* 

156. 


2. Sample 

Our s ample is based on updat ed version of the supplemental ta¬ 
ble of ICarilli & WalteJ (l2()13b . which contains all published CO 
line measurements at z > 2 by December 31, 2012. The updated 
version also include some submm sources published in 2013. We 
created a large sam ple of high-z submm g alaxies by adding 32 
sources observed bv iBothwell et al.l (l2013l) . In total, our sample 
contains 195 galaxies, detected at least in one CO transition, with 
a measured FWHM. 

3. Analysis 

We compute using the standard definition 

L'co = 3.25 X 10'^Sco^vv^Dld + (1) 

where S co^t' is the velocity integrated line flux, is the line 

luminosity in K Km s“' pc^, Vots is the observed central fre¬ 
quency of the line and is the luminosity distance in Mpc. 

If the molecular gas emission comes from thermalized, op¬ 
tically thick regions, is constant for all J levels. However, 
in reality, the ratio of L’^gd - 0) to higher order J - d 
depends on temperature and density. We use a verage ratio of 
J - dlL'ggd - 0) provided in Table 4 of iBothwell et al.l 
( |2()13|) . to calculate L'^gd - 0) from higher order CO observa¬ 
tions. Higher redshift galaxies are often observed only at mid-to- 
high-J CO transitions. We discuss uncertainties associated with 
this conversion in Section|4]in detail. 

It is always a concern with high-z objects, that they might be 
gravitationally magnified. Submm galaxies are not an exception. 
L'gg thus also needs to be corrected for a magnification factor 
(p), when the galaxy is magnified by a gravitational lens. We col¬ 
lected yu-values from the literature whenever available. However, 
as fi must be estimated from detailed modeling of the lens po¬ 
tential, a reliable measurement is not always available. We will 
examine how the correlation depends on the availability and am¬ 
plitude of ji in Section |4] but in brief, they do not change our 
main results. Most likely, majority of galaxies with missing p are 
not strongly magnified; otherwise a lens galaxy would have been 
detected and p would have been estimated. We present L'ggjp, 
whenever p is available hereafter. 


4. Results 

4.1. Lgg/p-FWHM correlation 

In Fig.l, we plot CO line width FWHM against L'ggjp for all 
galaxies in our sample. The solid line is the best-fit power-law, 
which is 

\og{FWHM) = (0.24 + 0.03) x \og{L'gglp) + (0.12 + 0.28) (2) 

The significance of the non-zero slope is 8.7 cr. Spearman’s cor¬ 
relation coefficient is 0.55 with p-value of 1.7x10 Kendall’s 
correlation coefficient is similarly significant as shown in Table 
1. The blue dashed line connects the median points in each bin. 
Errors on each point are small due to the large number of galax¬ 
ies in our sample. The median values increase with increasing 
L'^glp. These results show there exists a weak, but very signifi¬ 
cant correlation between L'ggjp and FWHM. 

The presence of such a correlation is not surprising. L'^glp 
is sensitive to the total mass of molecular gas, while the FWHM 
is sensitive to dynamical mass. The relation is similar to Tully- 
Fisher relation, which is also a power-law, in the local Universe 
(iTullv & Fisheij[l97^ . A large scatter around the relation is also 
expected due to the inclination effect, and the unknown p for a 
subset of galaxies. We measured a dispersion around the best- 
fit line of 43%, which is consistent with what one expects from 
randomly oriented thin disks (45%), suggesting p is close to 1 
for most of the galaxies with unknown p. 

4.2. Uncertainty in converting different transitions 

We test if conversions from different CO transitions can mimic 
the correlation. In Fig.l,we used different symbols for different 
CO transitions (4-3,3-2,2-l, and 1-0). In Table 1, we show cor¬ 
relation coefficients for these. For CO(4-3), we do not see a pos¬ 
itive correlation, while we see weak correlations for the other 
transitions. It is reassuring to see the correlations within single 
CO transitions, because they cannot be due to errors in the tran¬ 
sition conversion. The reason why we do not see a positive corre¬ 
lation for the CO(4-3) transition is not clear, but we point out that 
this is the smallest subsample, and that since a single transition 
only spans a small redshift range (and thus, a small luminosity 
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Fig. 2. is plotted against FWHM. We separated the sample into 

z>3 (red triangles) and z<3 (blue diamonds). The best-fit power laws 
for the two samples are; FWFIM (km/s) (z<3, blue) , 

and FWHM (km/s) oc (z>3, red). The slopes are consistent 

with each other, and with the slope derived from the full sample in Fig. 1. 
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Fig. 3. Luminosity distance in Mpc against z. The red circles and 
black triangles are our CO samples separated at z=3. The black dots 
are 192 supemovae used to calibrate our low-z sample. The black solid, 
green dotted and blue-dashed lines are (Dm, Oa)=( 0.3,0.7), (0.3,0.0), 
and (1.0,0.0), respectively. 

range), it is difficult to examine a correlation with a small sample 
in a single transition. It will be necessary to repeat this test using 
the wider frequency coverage of ALMA. On the other hand, we 
see signihcant correlations in all cases when we exclude a cer¬ 
tain CO transition from the sample. These tests suggest that the 
LJ-q-FWHM correlation is not an artifact of CO order conver¬ 
sions. 

4.3. Dependence on magnification 

Among our sample, 85 galaxies have known p, while the other 
110 galaxies do not. We should hrst note that many submm 
galaxies with unknown p are probably not strongly magnihed; 
otherwise we would have seen the lens galaxy/cluster and mea¬ 
sured p. With this in mind, we still test if an unknown p can 
mimic the correlation. In FiglT] we mark two extreme samples of 


galaxies; galaxies with ju > 3 in cyan, and those with unknown 
p in green. There is no apparent deviation from the best-ht re¬ 
lation, except that more ps, are unknown at higher redshift. In 
Table 1, we separate the sample into those with known p, and 
those with unknown p, both of which show a signihcant corre¬ 
lation. Because brighter submm galaxies are more likely to be 
lensed (due to the steep luminosity function), the correction to 
unknown p will steepen the correlation, i.e., the correlation will 
not disappear with the correction. It is re-assuring to see a posi¬ 
tive correlation using many submm galaxies with unknown p. 

4.4. Dependence on gaiaxy types 

We also separate the sample according to galaxy types in Table[T] 
and FigH] Unfortunately, 24pm sources and BzK galaxy samples 
are too small to examine the correlation. Both submm galaxies 
(SMG) and QSOs have a signihcant correlation, with submm 
galaxies having a greater signihcance. 

4.5. Redshift dependence 

Next we examine if the L^q-FWHM correlation changes with 
redshift. Unfortunately, due to the Hux limits of the telescopes, 
the higher redshift sample includes more luminous sources, 
while lower redshift sources have smaller luminosity due to vol¬ 
ume effects. With this limitation in mind, we separated the sam¬ 
ple into high (z>3) and low (z<3) redshift samples in Fig.2. The 
best-ht lines to each sample agree with each other within the er¬ 
rors, consistent with no signihcant redshift-evolution. The slopes 
are also consistent with each other, and with that of the full sam¬ 
ple in Fig.l. 

A concern for low level CO transition is that CMB tempera¬ 
ture becomes higher at the highest redshift end. For example, the 
CMB temperature at z=6 is ~19K, which is higher than the tem¬ 
perature needed to excite the CO(2-l) line (16.6K). Fortunately 
at high redshift, ALMA can only observe higher excitation lines, 
which are less affected. At z > 6, only CO(6-5) (7'(.=116.16K) 
and higher excitations can be observed with ALMA. The relation 
would be still valid at lower redshift (z ^ 5). However, the effect 
of the higher CMB temperature need to be examined carefully at 
higher redshifts. 

Although we did not hnd any sign of redshift evolution, we 
also note that it is expected the gas fraction in a galaxy should 
increase as a function of redshift. Therefore, it is important to re¬ 
peat this test when larger sample become available in the future. 

4.6. Estimating gravitationai magnification using 
L'i^qI p-FWHM reiation 

If we assume the LJ^^/yU-FWHM relation, the equation can be 
used to estimate p. We re-write the equation as follows. 


p^ (3) 

The units for and FWHM are K Km s ' pc^ and Km 
s ', respectively. The equation can be useful when, for example, 
a submm galaxy is detected behind a galaxy cluster, and subject 
to gravitational magnihcation, but the detailed lens model is not 
available. Using the equation|3] one can estimate p and thereby 
intrinsic luminosities of the submm galaxy. 
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5. Discussion 
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5.1. A possible cosmological application 

If the correlation is independent of redshift, the luminosity dis¬ 
tance to CO emitting galaxies can in principle be measured by 
estimating from FWHM. Such a distance measure would be 
useful because alre ady now, galaxies at z~6 have been de tected 
in CO with ALMA dRiechers et al.l20r^lWang et al.l201^ . This 
method can be used up to such a high redshift, where it is diffi¬ 
cult for other distance measures to reach. Here we assess whether 
such an approach could be feasible. 

in Figs 1 and 2 assumes a cosmology; 
(Qm, Oa, h)-(0.3,0.1,0.7). Therefore, we need to calibrate 
the correlation by using CO detected galaxies with a measured 
distance from another method. This, however, is difficult 
because CO detected galaxies are mostly at very high redshifts, 
where other distance measures cannot easily reach. 

There are a number of ways to calibrate the distances, 
but here we use a method often used to calibrate gamma ray 
burst. We use an empirical formula for the luminosity distance 
as a function of redshift based on 192 type la supernovae at 
0.16 8<z<1.755 (DillO^' ^cm = 6.96zi™ -t 14.79z' °2; equation 
1 of iTsutsui et al.l (l2009l) . This empirical formula does not as¬ 
sume any cosmology, but just that the type la supernovae are 
standard candles. We use the formula to assign luminosity dis¬ 
tances to our submm galaxies at z<3, and thus to calibrate our 
Lpg-FWHM correlation. Then, we apply the L'^g -FWHM cor¬ 
relation to measure luminosity distances to the z > 3 submm 
galaxies. 

The results are shown in Fig.3. Along with the data 
points, we plot three different cosmological models with 
(Qm, f^A)=(0.3,0.7), (0.3,0.0), and (1.0,0.0) in the black solid, 
green dot ted and blue-dashed lines, respectively. At this stage, 
it is clear that the scatter is more than an order of magnitude 
larger than needed to obtain meaningful constraints on the cos¬ 
mological parameters. 


5.2. Sources of scatter in the relation, and prospects for 
reducing it 


The scatter around the L^^-FWHM correlation (43%) is con¬ 
sistent with what is expected from randomly oriented discs, so 
if most submm galaxies are rotating spiral-like galaxies, then 
the largest uncertainty likely comes from the unknown inclina¬ 
tions. With the full operational array, ALMA’s spatial resolu¬ 
tion will be as good as 0.037” at 110 GHz (for CO observation) 
in the extended configuration. At z~6, this corresponds to 0.2 
kpc, which is sufficient to re solve submm galaxies, which typi¬ 
cally extend over several kpc (IWang et alj2013l : lDe Breuck et alJ 
l2014li . With such observations it will also be possible to exclude 
mergers and pressure-supported galaxies which may be causing 
some of the extreme outliers from the relation, and to select sam¬ 
ples of rotationally dominated galaxies. Furthermore, because 
the inclination is a random effect and its distribution is known 
(with <sin \>~0.19 and Icr ~0.22), one can statistically correct 
for the effect. For the Tully-Fisher relation, for example, it has 
been shown that even when no information on the inclination 
is available, by using maximum likelihood estimation, the true 
relation can be recovered with only 1.5 times larger statistical 
error than when inclinations are known with zero uncertainty 
(lObreschkow & Mevei 2013h. As a reference, using the Tully- 
Fisher relation, Sakai et al.l (l2000h obtained Hq- 1 l±A (random) 


+7 (systematic), i.e., 6 and 10% of statistical and systematic er¬ 
rors, respectively. 

Another possible source of scatter is calibration errors, as 
even state of the art sub-mm telescopes can only be calibrated to 
10% at best for individual observations. Such errors are also ex¬ 
pected to improve drastically as large, homogeneous statistical 
samples become available from ALMA, as the error on the me¬ 
dian reduces with Viv. In a sample of 10000 submm galaxies for 
example, the statistical error on the median will reduce to 1%. 

Note that this assumes the calibration error is random. If 
there is a systematic calibration error we are unaware of, this 
error will be added to the distance measurement. 

Another possible source of scatter in the relation which will 
be brought down in future ALMA surveys, is conversions be¬ 
tween CO transitions. In such samples it will be possible to 
study galaxies in the same transition as a function of redshift, 
or subsets of galaxies with well-constrained transition conver¬ 
sions. Once the observational and systematic errors have been 
brought down sufficiently to measure the intrinsic scatter of the 
relation it will to asses its possible cosmological applications. At 
the moment, it is impossible to tell. 
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